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ABSTRACT 

Massive star clusters observed in galaxy mergers are often suggested to be progenitors 
of globular clusters. To study this hypothesis, we performed the highest resolution 
simulation of a gas-rich galaxy merger so far. The formation of massive star clusters 
of 10 5 to 10 7 Mq, triggered by the galaxy interaction, is directly resolved in this model. 
We show that these clusters are tightly bound structures with little net rotation, due 
to evolve into compact long-lived stellar systems. Massive clusters formed in galaxy 
mergers are thus robust candidates for progenitors of long-lived globular clusters. The 
simulated cluster mass spectrum is consistent with theory and observations. Tidal 
dwarf galaxies of 10 8-9 Mq can form at the same time, and appear to be part of a 
different class of objects, being more extended and rotating. 
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1 INTRODUCTION 

Globular clusters (GCs) are an important fos sil record of 
the e volution of physical conditions in galaxies (Bek ki et al.l 
2008), because the formation of massive clusters is triggered 
by shocks and high pressures in the inter stellar medium 
; |van den Berghlll979l ; fAshman fc Zepj|200ll ). Fundamental 
properties along the Hubble sequence are a higher frequ ency 
of GCs around elliptical than disk galaxies (Harris 199j]) and 
a bimodal population in particular around early-type galax- 
ies, with a population of low-metallicity GCs and a popu- 
lation of younger, higher-metallicity ones (|Ashman fc ZepJ 
1 19921 ). 

Since the work s of ISchweizerl (|l987l ) and 
IWhitmore et all (|l993h . there has been increasing evi- 
dence that young massive star clusters (YMC) form in 
interacting and merging galaxies, and these are often 
proposed to be GC progenitors. Theoretically, massive 
clusters form in mergers because of shocks and high 
turbulent pressure in interacting galaxies, which favors the 
formation of tightly bound clusters rath er than unbound 
associations (|Elmegreen fc Efremovlll997l ). If GCs can form 
this way, given that most elliptical galaxies formed by 
galaxy mergers (Naab & Burkert 2003, Bournaud et al. 
2005), this mechanism could account for the youngest and 
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most met allic GCs that are mor e frequent around early- type 
galaxies (|Ashman fc Zepj|200ll ). These GCs would come in 
addition to those formed early in the U niverse as a result o f 
thermal instabilities in proto-galaxies (IFall fc Reesl 1985). 
strong shocks at the epoch of the Reionizat ion llCenlj2 001). 
or stripping of nucleated dwarf galaxies (|Freemanl Il990l ; 
ICao et aTll2007l ). Nevertheless, that an important popula- 
tion of GCs formed i n galaxy mergers is still challenged 
by some observations (jSpitler et al.ll2008l ). Whether or not 
YMCs formed during mergers contribute to the present-day 
GC populations is still an open question also bec ause the 
long- term evolution of such YMCs is uncertain fsee lde Griisl 
120071 ). A requirement for YMCs to become long-lived GCs 
is that they should be gravitationally bound, which is 
difficult to assess observationally, as the velocity dispersion 
of a cluster does not directly trac e its mass because of 
mass segreg ation l|Fleck et al.l l200rj ). The survival issue is 
further complicated by the the mass-loss of clusters if their 
IMF is too shallow, and by the tidal field of the parent 
galaxy which may progressively disrupt these clusters 
jMiocchi et~aill2006l ). 

Numerical simulations are a powerful tool to study 
galaxy mergers, but resolving the formation of star clus- 
ters in self-consistent models requires a huge dynamical 
range, because structures smaller than 100 pc should be re- 
sol ved in a simulated volu me larger than 100 kpc. M odels 
bv lBekki etakl (|2002D and iKravtsov fc Gnedinl <|2005l ) have 
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shown that the pressure and density required to form mas- 
sive bound star clusters could be reached in galaxy mergers. 
They could ho wever not identify individually forming GCs. 
iLi et al l (|2004h studied individual GC formation in an indi- 
rect manner, using a model where absorbing sink particles 
are assumed to form above a chosen density threshold. This 
allows the mass and spatial distribution of the putative GCs 
to be studied. However, the absorbing nature of sink parti- 
cles in such a model dooms any stellar association to be 
endlessly bound, so that whether or not long-lived bound 
objects can actually form in mergers is not probed by such 
a model. 

In this Letter, we present a simulation of a wet galaxy 
merger, which to our knowledge is the highest resolution 
model of this kind so far. We directly resolve the formation of 
dense structures resembling Super Star Clusters (SSCs) with 
typical masses of 10 5-7 Mq. Comparing their gravitational 
and kinetic energy, we argue that they are tightly bound 
and likely progenitors of long-lived GCs. 



2 NUMERICAL SIMULATIONS 

We use a particle-mesh code to compute the gravitational 
dyna mics of gas, stars, and dark matter l|Bournaud et al.l 
2005). The grid cell size is 32 pc up to 25 kpc from each 
galaxy center, 64 pc up to radii of 50 kpc, and 128 pc 
at larger radii. The density is computed with a Cloud-in- 
Cell interpolation, and a FFT technique is used to com- 
pute the gravitational potential, with a Plummer softening 
length of 32 pc for all types of particles. Gas dynamics is 
modeled with a sticky-particle scheme with elasticity pa- 
rameters /3t=/3 r =0.6. The star formation rate is computed 
using a Schmidt-Kennicutt law: the star formation rate is 
proportional to the gas density in each cell to the expo- 
nent 1.5. Gas particles are converted to star particles with 
a corresponding rate in each cell. Energy feedback from 
supernovae is accounted for with the scheme proposed by 
iMihos fc Hernquistl (|l994h . Each stellar particle formed has 
a number of supernovae computed from the fraction of stars 
above 8 Mq in a Miller-Scalo IMF. A fraction e of the 
10 51 erg energy of each supernova is released in the form 
of radial velocity kicks applied to gas particles within the 
closes t cells. We use e = 2 x 10~ 4 , as IMihos fc Hernquistl 
suggest that realistic values lie around 10 4 and less 
than 10" 3 . 

We use a total number of particles of 36 millions, cor- 
responding to 6 million particles per galaxy per compo- 
nent (gas, stars, and dark matter). The particle mass is 
3.4 x 10 4 Mq for initially present stars and 7 x 10 3 Mq for 
gas and stars formed during the simulation. To our knowl- 
edge, this is the highest-r esolution simulat i on of a wet galaxy 
merger published so far. Wctzstci n et alj l|2007h performed 
a merger simulation with 4 x 10 6 particles per galaxy, but 
only 45,000 particles for the gas comp onent. The s imula - 
tion of early-type galaxy formation by iNaab et al.l (|2007l ) 
uses 8 x 10 6 particles, but the large volume modeled in this 
cosmological ru n limits the resolution (softening) to 250 pc. 
I Li et al.l (|2004l ) had 5 x 10 5 gas particles per galaxy, which 
is one order of magnitude below our number and gives an 
equivalent particle mass only if one assumes a much lower 
galactic mass; they use a softening of 10 pc for gas particles 



but 100 pc for stellar particles, so that gravity is only mod- 
eled accurately at scales larger than 100 pc since the gravity 
of gas and stars cannot be decoupled (e.g.. |joslll996l ). 

The initial setup was chosen to be representative of an 
equal-mass wet merger at low or moderate redshift (O^z^l) 
with a gas fraction of 17%. Each spiral galaxy is made- 
up of a stellar disk with a Toomre profile of scale-length 
4 kpc truncated at 10 kpc, and a gas disk of scale-length 
8 kpc truncated at 20 kpc, which models the extended HI 
disks detected well beyon d the optical edge in most spirals 
l|Roberts fc Havnesl fl99i ) . The bulge has a Plummer pro- 
file of scale-length 0.8 kpc. The dark halo has a Burkert 
profile with a 8 kpc core radius, truncated at 75 kpc. The 
initial stellar mass is 2 x 10 11 Mq, 18% of this mass being 
in the bulge. The gas disk mass is 4 x 10 10 Mq. The dark 
halo mass is 5.5 x 10 11 Mq. The orbit in the merger model 
was chosen to avoid peculiar situations like coplanar disks 
or head-on encounters. It is prograde for one galaxy, with 
an angle between the disk and orbital plane of 30 degrees. 
It is retrograde for the other galaxy, and the angle between 
its disk and the orbital plane of 70 degrees. The pericenter 
distance is 25 kpc and the velocity at infinite distance is 
150 km s" 1 . 



3 RESULTS 

3.1 SSCs in isolated and merger models 

The evolution of the merging galaxy pair is shown on Fig. [T] 
Both galaxies develop tidal tails shortly after the pericen- 
ter passage while their central bodies merge together in less 
than 500 Myrs. One billion year after the first pericenter 
passage of the two spiral progenitors, the central elliptical 
galaxy is already relaxed, while tidal debris are still visible 
around it. As in all mergers of this type, the remnant has the 
shape, profile and kinematics of an elliptical galaxy. At the 
end of the simulation, the morphology resembles that of ob- 
served late-stage mergers, like IC 1182 on the last snapshot 
of panel (a) or NGC 7252 on panel (b) in Fig. [1] 

Shells and filaments are present as expected for a young 
merger remnant. Several compact and dense structures are 
visible in and around the elliptical remnant, over a large 
radial range. The most massive are manifest in the total 
stellar mass distribution, and the others are emphasized by 
plotting the density of the "young" stars alone (defined as 
those formed after the first pericenter passage). More than 
one hundred objects of this type are visible, with stellar 
masses m the 10 5 " 7 M range, diameters between 10 and 
100 pc. They are made-up almost exclusively of "young" 
stars and as shown below are gravitationally bound: they 
should therefore be considered as SSCs formed during the 
galaxy merger. We will show later that the three most mas- 
sive objects, above 10 s Mq, should be considered as systems 
of another class, namely Tidal Dwarf Galaxies. 

These SSCs are, for a large part, well above the mass 
resolution of our models: the number of particles describ- 
ing a cluster of 10 6 Mq is 140; the mass spectrum of the 
clusters formed during the merger is shown on Fig. [2] Thus, 
except for those of ~ 10 5 Mq, these SSCs are not rem- 
nant objects that were prevented to fragment because of 
a too small number of particles. To prove that the forma- 
tion of these SSCs is really triggered by the galaxy merger, 
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Figure 1. (a) Sequence of four snapshots showing the merger evolution. Time is indicated in Myr after the first pericenter passage. The 
total stellar mass density is shown - (b) Surface density of "young" stars (denned as those formed after the first pericenter passage) at 
the end of the wet merger simulation, 950 Myr after the first pericenter passage. Compact SSCs are seen all around the central elliptical 
galaxy, together with three massive tidal dwarfs in the remnants of the two long tidal tails. 



and not simply by a general instability in our galaxy mod- 
els, we show on Fig. [3] the isolated disk evolution at the 
same resolution as our merger model. It has been evolved for 
the same duration as the merger and the figure shows the 
"young" stars, defined as the stars formed after t = 300 Myr 
(the pericenter in the merger model). This model does not 
show many dense and compact substructures. A few asso- 
ciations of young stars are visible in the outer spiral arms, 
but none reaches a mass of 10 6 M and they are not as 
compact as the structures formed during the merger; these 
are only short-lived associations of young stars. This shows 
that our isolated disk galaxy model is overall stable against 
the formation of SSCs. This is consistent with observations 
of spiral galaxies where you ng SSCs form mostly in the in- 
ner resonant rings of bars l|Galliano et al.| [2005), while no 
such ring is formed in our isolated galaxy model over the 
Gyr-long evolution period. The formation of dense, SSC-like 
objects in our merger model is a direct consequence of the 
galaxy interaction itself. To check that the results are not 
crucially dependent on the sticky particle parameters, we 
ran another simulation with / 9 t =/3 r =— 0.6 instead of +0.6: 
colliding clouds pass through each other instead of bouncing 
back, which is a more dissipative case. The final system has 
similar SSC properties (Figs. [5] and [3}. 

A noticeable property of the interstellar gas during the 



galaxy interaction is its increased velocity dispersions. At 
t=450, which corresponds to the star formation peak, the 
average dispersion of gas clouds is 16 km s _1 and 17% of gas 
clouds experience local dispersions larger than 30 km s _1 . 
In the isolated disk at same instant, the average value is 
11 km s _1 , and only 3% of the clouds are above 30 km s _1 . 
This is consistent with the large dispersio ns observed in in- 
teracting systems l|Elmegreen et al1ll995l ). As a result, the 
Jeans mass is increased and more massive bound structures 
can form in the gas. Nevertheless, such massive structures 
could become unbound after star formation has operated 
(|Gever fe BurkertJl200 if ). Actually, the gas pressure is en- 
hanced during the galaxy merger, with local spikes of high 
pressure in the caustics and shocks caused by the interac- 
tion: the average pressure is increased by 55% and the 95th 
percentile of the pressure distribution by a factor 9. Star for- 
mation in such high-pressure regions has a high efficiency, 
so that the resulting stellar systems can remain bound: this 
enables the formation of massive bound star clusters. 



3.2 Evolution into globular clusters 

The young SSCs observed in galaxy mergers and reproduced 
in our model are often suggested to be progenitors of GCs. 
The infant mortality can affect these SSCs but a large frac- 
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Figure 2. Mass spectrum of the stellar structures formed in the 
merger model. The dashed histogram is for the test simulation 
with reversed sticky particle parameters. Power-law spectra with 
slopes -2 (dashed) and -3 (dotted) are shown. 



10 kpc 





Figure 5. Stellar velocity and dispersion fields for 2 TDGs and 
3 SSCs. The stellar mass and typical rotation velocity to disper- 
sion ratio (V/cr) are indicated. TDGs are rotating disks with V/cr 
larger than 1. GC progenitors are dominated by random motions. 



Figure 3. Isolated galaxy evolved over the same duration as the 
merger. The young stellar density is shown and does not harbor 
substructures as dense and compact as in the merger remnant. 

tion can still be long-lived (|de Griisll2007M and may thus be- 
come GCs. This however requires that they are dense enough 
and tightly bound by their own gravity, which our model al- 
lows to check directly. Measuring the half-mass radius of 
the SSCs in our model, we find that more than two thirds 
have radii between 10 and 35 pc: even if they are above the 
mass resolution limit, they are around the spatial resolu- 
tion, which suggests that they have reached an equilibrium 
because of the softening of Gravity at short distances. We 
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Figure 4. Distribution of l/2Ea + Efc for the SSCs in our model, 
where Eq is computed for a true, non-softened Gravity. This value 
would be (— oo in logscale) for virialized systems. Our SSCs 
are characterized by Eq/2Eh < —1, which indicates that their 
collapse was stopped by the spatial resolution limit; they would 
become even more compact at infinite resolution. 



have measured for each cluster the kinetic energy Ek and 
gravitational energy Eg for the real, non-softened gravity. 
While the value of 1/2Eq+Ek should be zero for structures 
in pure gravitational equilibrium, we find that this quantity 
is significantly negative (equivalently Eq/2Ek < — 1). This 
confirms that the clusters have reached a pseudo-equilibrium 
resulting from the softening, and would become more com- 
pact at infinite resolution, with half-mass radii then signifi- 
cantly smaller than the 10-35 pc found here. 

The SSCs formed during the galaxy merger are thus ro- 
bust candidates for GC progenitors. They have masses typi- 
cal of GCs down to our mass limit. They are compact objects 
with typical half-mass radii that would be smaller than 10 pc 
at infinite resolution. They are bound and long-lived: many of 
them are already formed 500 Myr before the final snapshot 
and their boundedness would increase at higher resolution. 

3.3 Mass spectrum and tidal dwarf galaxies 

The mass spectrum of young stellar objects is shown on 
Fig. [2] Masses of 10 Mq can be affected by the mass res- 
olution of our models, and the turnoff of the mass func- 
tion below 5 x 10 Mq cannot be considered as robust. At 
higher masses, the mass spectrum decreases with a power- 
law shape. 61 objects have masses of 10 5-6 , 66 of 10 6-7 , only 
6 of 10 7-8 . There is a noticeable gap between the most mas- 
sive SSCs and the 3 objects with masses of 10 8 " 9 M n , which 
resem ble Tidal Dwarf Galaxies (TDGs, iBournaud fc Duel 
2006). TDGs are much less concentrated than the SSCs: 
SSCs of a few 10 6 M have half- mass radii around 20 pc, the 
nine most massive SSCs of 1 to 5 x 10 7 Mq have an average 
radius of 48 pc, while the three TDGs have radii of 1-2 kpc. 
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This leads to typical surface densities (oc M/R 2 ) one or- 
der of magnitude lower for the TDGs. While the SSCs/GC- 
progenitors are compact structures supported by random 
motions with little net rotation and V/a ratios much smaller 
than 1, TDGs have a large angular momentum (V / /cr~2-3) 
and take the form of rotat ing gas disks (see Fig . [5} as also 
found in observations (e.g. iBournaud et al.ll2007l ). 

Because of these physical differences, in addition to the 
possible gap in the mass function, our model suggests that 
two families of objects can form in galaxy mergers. Compact 
SSCs, in which rapid star formation occurs, are robust can- 
didates for GC progenitors. TDGs are not simply the high- 
mass end of the SSC distribution: they are made-up of ro- 
tating disks forming stars at a more mode rate pace, and will 
likely evolve into dwarf sate llite galaxies (jBournaud fc Dug 
120061 : iMetz fc Kroupal [2007b . The SSCs/GCs form locally 
by the gravitational collapse of overdense regions in the tur- 
bulent interstellar medium, following the mechanism stud- 
ied for instance by IWetzstein et al.l ifeoOTl ). while TDGs arc 
formed in regions where large amounts have been accumu- 
lated during the interaction o r in t he p rogenitor d i sks, a s 
proposed bv lElmegreen et all (|l993h and iDuc et al.1 l|2004h . 

As for the mass spectrum of SSCs themselves, above 
the mass limit and up to a few 10 7 Mq, the best-fitting 
slope in a iV cx M a model would be a ~ —3, as shown 
on Fig. \%\ However, the -2 slope exp ected from theoretical 
models (jElmegreen fc Efremovlll997l ) is compatible with the 
bulk of the mass distribution, assuming that the spectrum 
is truncated around 10 7 Mq, which may indee d be realistic 
and c ould result from various physical effects (|Gieles et al.l 
2006). Models spanning a larger mass range or repeated runs 
to increase the statistics are be required to firmly assess the 
shape of the cluster mass function. Nevertheless, our present 
model is compatible with observations s uggesting a slope 
around -2 or steeper (|Weidner et al.l 



4 SUMMARY 

Using a high-resolution simulation of a galaxy merger, we re- 
solve the formation of structures down to masses of 10 5 M©. 
This enables us to directly reproduce the formation of nu- 
merous Super Star Clusters, as observed in merging galaxies. 
We find that these are dense, tightly bound structures that 
will evolve into compact stellar systems and are likely pro- 
genitors of numerous globular clusters (GC). Tidal dwarf 
galaxies of a few 10 8 Mq are also formed in our model; 
they appear to be a different type of objects, less concen- 
trated and resolved as rotating disks. Our results suggest 
that both TDGs and SSCs can form at the same time in 
galaxy mergers. The mass function of GC progenitors in 
this model has a power-law shape with a slope between -2 
and -3, likely truncated around 10 7 Mq. The formation of 
globular cluster is much more efficient in merging galaxies 
than in isolated disk models. The total mass of the GC pro- 
genitors in our merger model amounts to 4% of the available 
gas mass, which is 0.7% of the total baryonic mass. The ef- 
ficiency could be even higher with the large gas fractions 
observed in high-redshift disks (|Daddi et al.ll200"8T ). GC pro- 
genitors form in the tails and caustics with high turbulent 
speed and pressure around the central merger remnant. The 
material in such regions has been expulsed from the spiral 



disks where it has been previously enriched; it is known t o 
have relatively high metallicities (Wcilbac her et ail 120031 ). 
so the GCs formed this way should also be metal-rich. The 
excess of GCs around elliptical and lenticular galaxies com- 
pared to spirals could then be explained by this mechanism, 
given that mergers can form early-type galaxies and GCs at 
the same time. 
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